Moil. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN style file v2.2) 



On the association of G343. 1-2.3 and PSR 1706-44 

R. Dodson \ K. Golap ^ 

^ School of Mathematics and Physics, University of Tasmania, GPO Box 252-21, Hobart, Tasmania 7001, Australia; 
Richard. Dodson@utas. edu. au 

^ National Radio Astronomy Observatory, P.O. Box 0, Socorro, NM 87801 
kgolap@nrao. edu 



1 February 2008 



ABSTRACT 

The association of G343.1-2.3 and PSR 1706-44 has been controversial from 
its first proposal. In this paper we address the difhculties, and argue that the 
association is still likely. New evidence comes from images of G343. 1-2.3 ob- 
tained using the Australia Telescope Compact Array (ATCA), and the pulsar 
obtained using the CHANDRA X-ray observatory. Mosaicing was required to 
cover the full extent of G343.1-2.3, and we present the polarisation images 
from this experiment. Also an X-ray pulsar wind nebula has been found in 
the archived CHANDRA observations, with the correct morphology to sup- 
port the association. The ATCA observations confirm the much larger extent 
of the SNR, which now encompasses the pulsar. The X-ray morphology points 
back toward the centre of the SNR, indicating the direction of the proper mo- 
tion, and that the PSR and SNR are associated. 
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1 INTRODUCTION 

The pulsar PSR B1706-44 is one of only seven that are known to emit gamma-rays in the 
GeV energy range ( [I'hompson et al. 1992| ) and one of only three that have been detected in 
the TeV range ( IKifune et al. 19951 ; IChadwick et al. 19981) . [Becker et al. (19951) reported that 



the power law spectrum from the compact source at the site of the pulsar was best explained 
by synchrotron origin, presumably a nebula around the pulsar. Finley et al. (199^ ) found 
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that PSR B1706-44 is surrounded by a compact X-ray nebula of radius ~ 0.3 pc using ASCA 
data. This is also true of the other two TeV emitters, the Crab and Vela pulsars. The TeV 
gamma-rays, which appear to be unpulsed, can be explained as resulting from the Inverse 
Compton interaction of the high energy electrons, which produce the X-ray synchrotron 
emission, with ambient photons ( [Aharonian et al. 1997[ ). 

The original discovery that the COS-B gamma-ray source was a pulsar was made by 
[Johnston et al. (1992| ). The spin down luminosity is high, 3.4 x 10'^^ ergs, as would be 
expected. [McAdam et al. (1993| ) published a map made by the MOST telescope at 843 MHz 
of the area around the gamma-ray source. It showed a semicircular arc of emission, which 
has subsequently been denoted as the supernova remnant (SNR) G343. 1-2.3, with the pulsar 
seemingly embedded in a small feature at its south eastern extremity. It was argued that 
the approximate distance of 3 kpc for the remnant, derived from the surface brightness- 
diameter relationship (S-D) was compatible with the pulsar's dispersion measure, although 
the distance indicated by the widely used interstellar electron density model of [Taylor ^ 
Cordes (1993|) would be only 1.8 kpc. The E-D relationship was based on the flux density 
values from single dish observations, as both the VLA and the MOST integrated flux fell 
short of what would be expected, implying that the broad structure was being resolved 
( IMcAdam et al. 1993| ). 

If the association of the pulsar and the supernova remnant is real, a transverse velocity 
of ~ 900 kms~^ is required for the pulsar to have moved from the approximate geometric 
centre of the SNR arc in the characteristic spin down time of 17 kyr. This is high, but not 
the highest measured value for a pulsar, nor as high as that implied for some seemingly 
well-established pulsar-SNR associations ( [Manchester et al. 1991| ; [I'rail fc Kulkarni 1991[ ). 

Frail et al. (1994) imaged the area around the pulsar at 20 cm and 90 cm with the 
VLA and cast doubt on its association with G343. 1-2.3. These arguments were based on 
the morphology, and the fact that the dispersion measure (DM) and S-D distance were 
discrepant. The latter was subsequently countered by a 21 cm hydrogen line absorption 
measurement of the pulsar by [Koribalski et al. (1995[) which gives a distance range of 2.4 to 
3.2 kpc. Most recently the region has been imaged in more detail by [Giacani et al. (2001 ) 
at 20, 6 and 3.6 cm and they concurred with [Frail et al. (1994 ) description of a nebula size 
of 3.5' X 2.5'. 



Nicastro et al. (1996) also argue against the association. They report a maximum value of 
only 27 kms~^ for the magnitude of the transverse velocity of the pulsar calculated from the 
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interstellar scintillation. However in their later paper ( [Johnston et al. 1998| ), this estimate 
was revised up to 100 kms~^. The possibility of a direct determination of the proper motion 
of PSR B1706-44 has been investigated, but no phase reference sufficiently strong for the 
current Australian VLBI network (the LBA) could be found. 

We have made a number of high sensitivity, low surface brightness images of the remnant 
to determine its full extent. Furthermore we have used CHANDRA archived data to search 
for a pulsar wind nebula (PWN) to provide evidence for the pulsar interaction with the 
interstellar medium (ISM). 



2 RESULTS FROM THE COMPACT ARRAY 

The Australia Telescope Compact Array has been used to re-observe G343. 1-2.3, with a 
nineteen pointing mosaic. Mosaicing was required both to cover the extent of the source and 
also to assist in the recovery of the flux on the shortest uv spacings, which were missed in the 
MOST and the VLA observations. Primary calibration was against PKS1934-638 which has 
an integrated flux of 16.4 Jy at 1.384 GHz. The mean observing epoch was September 1998. 
Deconvolution used the mosaic version of the standard miriad maximum entropy methods 
for polarised images; pmosmem. The RMS noise in the final total intensity image is 0.6 mJy. 

Figure shows the total intensity images of the whole remnant (Figure |^ a) and of 
the region around the PSR (Figure b). The pulsar region is overlaid with the polarisation 
vectors. The full image is very similar to those from MOST ( |McAdam et al. 1993|) , and indeed 



has very similar resolution (70" x 47") (the longer baselines, to the 6km antenna, were not 
used due to the uv poor coverage). We have included as a constraint in the deconvolution 
the observed total flux for the SNR, 30 Jy, observed with the University of Tasmania's 26m 
antenna. 

When the total flux from single dish measurements is included, the images obtained are 
in agreement with the other observations, i.e. the extent agrees with the very low frequency 
images from the VLA ( |l^rail et al. T994D and MRT ( podson 1997| ; |Golap 1998| ), and the 



spectral index has a simple monotonic spectrum of -0.5 ( podson et al. 1999[ ) 



Figure ^ shows the polarisation angle overlaid on the polarised intensity. The polarised 
fraction is typically about 20 percent with a randomly orientated polarisation angle which 
changes rapidly around the bright ring. Wisps of polarised emission can be seen to extend 
beyond the originally defined ring, and throughout the region with broadscale emission. This 
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confirms that the broadscale emission is associated with the synchrotron emission from the 
SNR. 

The pulsar is a point source with a flux density of 10 mJy, in line with that found in the 
scintillation fluctuation observations of [Johnston et al. (1998| ), and is 70 per cent linearly 
polarised, and 25 per cent circularly polarised. The mean pulsar position is 17:09:42.71 - 
44:29:08.1 with an RMS of 0.5 arc seconds between epochs, within the errors estimates of 



Wang et al. (2000| ), but is slightly south of the VLA observations; that is -44:29:06.6. 

It is instructive to compare the VLA images of [Frail et al. (1994| , figure 4) with the image 
in Figure ^ b). The pulsar lies on a spur of emission which is clearly a continuous feature 
covering more than 30'. Thus it is extremely unlikely to be a cometary tail from the pulsar. 
Furthermore the width of the spur, 4', is too broad to be due to ram pressure of a moving 
pulsar. It is intriguing that this is approximately the size that would be expected from a 
static pulsar wind, derived from the Sedov equations with a continuous input of energy, 
and with the E of PSR1706-44 ( [Rees fc Gunn 1974| ). Taking the expressions from [Frail fc 
Scharringhausen (1997[) , the equation of the static shock radius (Rs) is; 



Rs = (^)'^'^'^'PC ^ ^-QEli'tTuo'^'D^^l arc min (1) 

Where i?34 is the spin down energy in unit of 10^^ ergs s~^, no is the number density, 
^5 is the age in 10^ years and D^pc is the distance in kpc. Using the standard values gives 
1.7'no ^''^ using the lower limit on the distance. 

The formula for a ram pressure shock driven PWN (Rw) is; 

R^ = ( . ^ 0A75EliW^'u^,lD,^\ arc sec (2) 

or 3.7"no ^^^z^isQ. 1^150 is the velocity in units of 150 kms~^. 

It is clear that the radius of the feature seen in the radio band cannot be due to the ram 
pressure shock unless the density is extremely low (2 x 10"*^). It is consistent with the size 
expected for a static PWN, with a number density of the order suggested by the analysis 
of the i^OS'^Tdata, i.e. 0.09 cm'^ ( podson 1991 ; Podson et al. 1999[ ; [Becker et al. 19951) . 
These observations are contaminated by the nearby low mass X-ray binary, LMXB 1705-44, 
but a best fit of a thermal spectrum to the X-ray emission from the SNR gives a distance 
of 3.1 kpc, and an age of 8.9 kyrs. These values are derived using the X-ray luminosity 
to distance model of [Kassim et al. (1994[) . The local density based on the pulsar age is 
0.09 cm~^, whilst that from the Sedov model age is 0.02 cm~^. 
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3 OBSERVATIONS WITH THE CHANDRA X-RAY TELESCOPE 



We took 48 ksecs of HRC data and 17 ksecs of ACIS data from the CHANDRA X-ray 
observatory archive and imaged the region around the pulsar. These observations were made 
in August and February 2000, with roll angles 180° apart. A similar low surface brightness 
extended feature is clearly visible in both the HRC and the ACIS data, as first reported by 
Finleyet al. (199^ ). Their ASCA observations, complicated by contamination from the bright 



LMXB 1705-44, saw evidence for a 14-33" compact nebula in the radial surface brightness. 
The feature as seen by CHANDRA as it is visible in both data sets, can not be due to the 
space craft dither, nor due to a soft proton flare or any other random and unrelated event. 

We have found the peak energy of the pulsar emission to be around 1 keV from the 
^C/5'data, and have 'CLEANed' the HRC data (smoothed to 0.3") with the PSF for that 
energy. CLEAN is a standard deconvolution technique where a small fraction of the PSF is 
subtracted interatively from the image, until a conversion condition is met ( Hogbom 1974 ). 
Our convergence criterion was satisfied when the peak emission was more than half the 
PSF width from the pulsar position. Significant emission (10 percent of peak) remained at 
this point. It was perfectly feasible to clean this feature iteratively as well, but equally good 
results were obtained by fitting a Gaussian model around the pulsar position. Figure § shows 
the emission after the cleaning the point source with the PSF; its position is marked with 
a cross. This confirms the object seen by Pinley et al. (1998]) , but finds a two dimensional 
structure and a smaller extent which can be explained by the degree of contamination in 
the ASCA image. The object has a very compact nebula (but much greater than the PSF) 
around the pulsar which is best fitted by a Gaussian of about l.liO.l", and a more tenuous 
trail pointing nor- north west, and extending about 5" back from the pulsar. The angle of 
symmetry of this tail is —15°, toward the SNR centre. The brightness of the tenuous nebula 
is about 10 times greater than that of the point spread function of the CHANDRA mirrors; 
i.e. 1 percent of the pulsar peak, whereas the PSF would be 0.1 percent. The most likely, 
and most natural, physical model of such a structure is a bow shock, with the more tenuous 
tail left behind in the wake. 

The pulsar to bow shock distance at the closest approach (1.1"), at 2.4 kpc, is just 0.05 
pc. If the scintillation velocity limit found by Johnston et al. (1998 ) is correct then the 
density is required to be greater than 25 cm~'^, or, if the densities are those found from the 
i? 05"^ T observations (Podson 1997 ; Dodson et al. 1999 ; Becker et al. 199"5|) , the velocities 



6 Dodson & Golap 

are in the range 1000-2000 kms^^, greater than that required to move from the centre of 
the SNR to the current position. More typical densities imply reasonable velocities (500 
kms~^ implies a density of unity). These are obviously much higher that those seen by the 
scintillation velocity, but we stress that this is an indirect measurement with a large intrinsic 
scatter, and only a VLBI observation will settle the question, for which the LBA will need 
to be upgraded. 

From the HRC image we cleaned out 294 counts for the pulsar before the convergence 
criterion was reached. The integrated X-ray flux from the compact nebula is 390 counts and 
the tenuous tail is 120 counts. We fitted a number of models to the ACIS data in sherpa 
2.2. The best fit to the source spectrum was found with joint fitting of the two sources; 
the region around the pulsar as a black body plus a power law portion and the PWN as 
an independent power law. These were attenuated with the usual column density absorp- 
tion. The fact that the same spectrum could not be fitted to both features confirms that 
it is not the result of incorrectly removed spacecraft dither. The model fits for the black 
body can be used to calculate the diameter of the emission region from the normalisation 
constant {228. 5tlllRlmD^okpc)- Using the Hi dist ance limits gives the radius of the emit- 
ting region as 3.6l!f 3-4.8I1 g km. The unabsorbed X-ray flux gives the range of intrinsic 
luminosities to be 5.4-7.1 xlO^^ ergs s~^ from the black body (presumably the pulsar) and 
7.6-9.9 xlO'^^ ergs s^^ from the power law regions (the PWN). The fits are shown in figure 
^. The reduced Chi-squared is 0.6 and a Q-statistic of 99%, implying that the errors have 
been underestimated. The column densities agree with those from the iJO^^T observations 
( IBecker et al. 1995|) . 

The fitted model parameters are; 
nn/lO^' Tsij/keV reduced x 

5.3 ±0.8 0.15 ±0.02 1.7 ±0.5 1.5 ±0.2 31/53 
The integrated X-ray flux between 0.5 and 10 keV for each component in ergs cm~^s~^ 



is; 
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Power Law, 



pwn 



Received 



0.7 X 10-13 1.6 X 10- 



4.6 X 10- 



Unabsorbed 5.8 x lO^^^ 2.2 x 10^ 



5.9 X 10 



-13 



On the association of G343. 1-2.3 and PSR 1706-44 7 

4 CONCLUSION 

Initially the arguments against the association were based on Kaspi's criteria ([Kaspi 19961) . 
These can be summarised in order of ease of observability as; same position, same distance, 
same age, velocity reasonable and proper motion correct. It was felt that the pulsar DM 
distance and the SNR S-D distance did not match, and the pulsar was at the edge of the 
SNR and had too low a velocity to have travelled from the geometric centre of the SNR. 
Also there was no morphological evidence to imply a velocity direction, and the scintillation 
velocity was too low. Since then work ( |Koribalski et al. 1995| ) has shown that the DM distance 
was incorrect, and our ATCA observations have found that the extent of the SNR has been 
underestimated. The Chandra data shows a morphological signature best interpreted as the 
path of the pulsar. The proper motion of the pulsar is the remaining question, and work is 
under way to attempt to answer that. 

It is clear from both the low total flux values found by the MOST and VLA maps 
that there is a smooth broad component to this SNR. The ATCA observations confirm the 
extent of the broad emission, and that this is non-thermal in character. Taking this smooth 
component as the true limits of the SNR enlarges the extent and moves it southwards. This 
places the pulsar PSR1706-44 within the SNR shell, rather than on the rim, thus satisfying 
Kaspi's first criteria. 

The X-ray PWN observed by CHANDRA is consistent with the pulsar travelling from 
the general direction of the SNR centre, with either a very slow velocity (which would agree 
with the scintillation results) or a low local density (in agreement with the i^OS'/lT results). 

We expect a SNR to be associated with PSR1706-44 because of its youth. The association 
with G343. 1-2.3 was rejected in the past, but we have shown that the grounds for the 
rejection are not compelling. The question of the pulsar velocity could be resolved by finding 
a phase reference for VLBI observation, with either the SKA or an upgrade to the current 
LBA. We will be able to rule out the upper velocity limits with the ATCA within the next 
few years but will not be able to measure the lower limits for a few decades. 

Finally the result stresses the importance of improving the ROSAT observations of this 
SNR with a mosaic over the whole area of the radio emission using either CHANDRA or 
XMM. 
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Figure 1. Supernova Remnant G343.1-2.3 at 1384 MHz. a) the full remnant and b) the region around the pulsar, with 
polarisation angles overlaid; a 10 mjy polarisation scale bar is shown 
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Figure 3. PSR1706-44 PWN in X-rays, from the Chandra HRC. Scale is in counts per doctor pixel, the 5, 10 and 20 a contours 
are shown. The FWHM of the used HRC PSF is shown at the site of the pulsar. 




Figure 4. Spectral fit to the pulsar region (black crosses) and the PWN (red squares). 



